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PATENT 

ATTORNEY DOCKET NO. : 10085US01 

COLOR PROCESSING 
3ackaround of the invention 
5 The invention relates generally to color processing 

and, more particularly, to techniques for color 
characterization and transformation. 

Since zhe introduction of the CIE (Commission 
International de I'Eclairage) color measurement system in 
10 the early 1930' s, many different color spaces have been 
proposed for different applications. A color space, also 
referred to as a color "metric, " is essentially a coordinate 
system by which a color can be quantified. 

A color space can be used to characterize the color 
15 output of a color imaging system relative to other color 

imaging systems. By characterizing multiple color imaging 
systems, the color space facilitates using different imaging 
systems to produce matching colors. An "ideal" color space 
would calculate a color mapping between different color 
2 0 imaging systems that achieves an acceptable color match 
between the systems without subjective or empirical 
adjustment . 

Color spaces differ in the parameters expressed on 
their coordinate axes and the manner in which the parameters 

25 are calculated. CIE color spaces use CIE Standard Observer 
functions that are based on color matching functions and 
result in a unique set of tristimulus values XYZ for any 
color measured under specified conditions. The tristimulus 
values XYZ are calculated from the spectral output of either 

30 an additive or subtractive color system convoluted with the 
response function of either a 2 degree or 10 degree Standard 
Observer. In the case of reflective hard copy, the spectral 
reflectance curve is typically convoluted with a standard 
illuminant to estimate the expected spectral output of the 

35 reflective color. 



One CIE color space is the CIELAB color space. In 
this color space, L* represents lightness, a* represents 
redness -greenness, and b* represents yellowness-blueness . 
The CIELAB color space employs a modified von Kries 
5 chromatic adaptation algorithm. According to the modified 
von Kries chromatic-adaptation transform, a description of 
which can be found in Gunter Wyszecki and W.S. Stiles, Color 
Science: Concepts and Methods, Quantitative Data and 
Formulae, section 5.12, John Wiley & Sons, Inc., 1982, the 

10 L*a*b* color spaces make use of white reference tristimulus 
data. The modified von Kries chromatic-adaptation transform 
involves dividing the tristimulus values XYZ obtained for a 
color produced by a particular color imaging system by white 
reference tristimulus values for the system. For example, 

15 the X, Y, and Z tristimulus values of the color under study 
can be divided, respectively, by the X, Y, and Z tristimulus 
values for a perfectly diffuse white reflector. Thus the 
von Kries approach defines both neutral and chromatic colors 
relative to the "white reference" representing the XYZ 

20 tristimulus values of the perfectly diffuse white reflector. 
The equations for the CIE 1976 CIELAB color space are as 



follows : 

L* = 116 x f (Y/YJ - 16 [1] 

a* = 500 x [f (X/X n ) -f (Y/Y n ) ] [2] 

25 b" = 200 x [f (Y/Y n ) -f (Z/Z n ) ] [3] 

f (w) = (oj) 1/3 cj>0. 008856 [4] 

f (w) =7.78 7 ((j) +16/116 wssO. 008856 [5] 



where X n/ Y n , and Z n are the tristimulus values of a 
perfectly diffuse white reflector under specified viewing 
3 0 conditions. The viewing conditions are determined by (1) 

the illuminant, e.g., D 50/ and (2) the Standard Observer (2° 
or 10°) . 
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Summary of the Invention 
In general, in one aspect, the invention provides a 
method of characterizing a color imaging system. The method 
comprises obtaining first data indicative of output of the 
5 color imaging system. The first data is processed, to yield 
second data, according to a color appearance model that 
varies in accordance with neutrality of colors indicated by 
the first data. 

In general, in another aspect, the invention 

10 provides a computer program product residing on a computer 
readable medium, for characterizing a color imaging system, 
and comprising instructions. The instructions are for 
causing a computer to obtain first data indicative of output 
of the color imaging system and to process the first data, 

15 to yield second data, according to a color appearance model 
that varies in accordance with neutrality of a color 
indicated by the first data. 

In general, in another aspect, the invention 
provides a method of producing a color on a device. The 

20 method comprises obtaining first data associated with a 

first device and indicative of a first color. Second data 
are determined that are related to stimulus data of the 
first device by a color appearance model that converts input 
data to output data using a white reference vector that 

25 varies in association with a neutrality of a color indicated 
by the input data. A second device is actuated according to 
the second data to produce a second color to approximate the 
first color. 

In general, in another aspect, the invention 

30 provides a computer program product residing on a computer 
readable medium, for producing a color on a device, and 
comprising instructions. The instructions are for causing a 
computer to obtain first data associated with a first device 



and indicative of a first color. The instructions are also 
for causing the computer to determine second data related to 
stimulus data of the first device by a color appearance 
model that converts input data to output data using a white 
5 reference vector that varies in association with a 

neutrality of a color indicated by the input data. The 
instructions are also for causing the computer to actuate a 
second device according to the second data to produce a 
second color to approximate the first color. 

10 In general, in another aspect, the invention 

provides a method of producing a color with an emissive 
device using absolute colorimetry. The method comprises 
obtaining first data indicative of a first color. Second 
data are determined that are related to the first data by a 

15 color appearance model that uses a white point of the 

emissive device as a white reference vector. The emissive 
device is actuated according to the second data to implement 
absolute colorimetry to produce a second color to 
approximate the first color. 

20 In general, in another aspect, the invention 

provides a computer program product residing on a computer 
readable medium, for producing a color with an emissive - 
device using absolute colorimetry, and comprising 
instructions. The instructions are for causing a computer 

25 to obtain first data indicative of a first color and to 

determine second data related to the first data by a color 
appearance model that uses a white point of the emissive 
device as a white reference vector. The instructions are 
also for causing the computer to actuate the emissive device 

3 0 according to the second data to implement absolute 

colorimetry to produce a second color to approximate the 
first color. 



- 4 - 



In general, in another aspect, the invention 
provides a method of characterizing an emissive device for 
absolute colorimetry. The method comprises obtaining first 
data indicative of output of the emissive device. The first 
5 data are converted to second data using a color appearance 
model that uses a white point of the emissive device as a 
reference white vector. The second data are provided for 
use in absolute colorimetric color reproduction. 

In general, in another aspect, the invention 

10 provides a computer program product residing on a computer 
readable medium, for characterizing an emissive device for 
absolute colorimetry, and comprising instructions. The 
instructions are for causing a computer to obtain first data 
indicative of output of the emissive device and to convert 

15 the first data to second data using a color space that uses 
a white point of the emissive device as a reference white 
vector. The instructions are also for causing the computer 
to provide the second data for use in absolute colorimetric 
color reproduction. 

20 In general, in another aspect, the invention 

provides a method of characterizing colors for reproduction 
between a first device and a second device. The method 
includes normalizing first tristimulus values indicative of 
a color of the first device using local black point values. 

25 The normalized first tristimulus values are transformed to 

obtain color values indicative of modified cone responses of 
the human eye. The color values are chromatically adapted 
from a local condition to a reference condition. The 
adapted color values are transformed to obtain second 

30 tristimulus values. 

Embodiments of the invention may provide one or more 
of the following advantages. The invention provides more 
accurate color characterization over the entire range of a 



color imaging system than color characterization using a 
fixed white reference point. The helps ensure substantial 
uniformity of the color characterization for both light 
color shades and more intense colors, and for neutral and 
non-neutral colors. Emissive devices can be caused to 
display colors,, both saturated and neutral, that match well 
to corresponding colors on absorptive devices. 

Other advantages will be apparent from the following 
description and from the claims. 

Brief Description of the Drawings 
FIG. 1 is a block diagram of a system, in accordance 
with the invention, for characterizing a color imaging 
system; 

FIG . 2 is a schematic diagram of the use of local 
and common white references depending on colors indicated by 
a color imaging system; 

FIG. 3 is a functional block diagram of the system 
shown in FIG. 1; 

FIG. 4 is a block diagram of a system, in accordance 
with the invention, for performing a color transformation; 
and 

FIG. 5 is a functional block diagram of the system 
shown in FIG. 4 . 

Description of Preferred Embodiments 
The invention provides methods and systems for 
characterizing color to accurately reproduce the color from 
one medium or color system to another where the two systems 
have different illuminant and/or media white points. 

As shown in FIG. 1, a system 10 for characterizing a 
color imaging system includes a processor 12, a color 
detector 14, a color data memory 16, a substrate 20, and a 
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L*a*b*. memory 28, Characterizing, also commonly referred to 
as "profiling, " can be used to analyze a single color 
imaging system, or as a basis to transform the color 
response of a "source" color imaging system to match the 
color response of a "target" color imaging system. The 
color response of a color imaging system is the correlation 
between a range of input color values (i.e., device stimulus 
data) and measured colors. The measured colors are produced 
according to the input values for color-producing systems 
such as displays or printers, and are the colors measured by 
color-measuring systems such as scanners. The device 
stimulus data can be, e.g., RGB or CMYK values. The 
stimulus data span ranges of discrete values of the system. 
For example, for an emissive RGB device, the values of R f G, 
and B can each range from 0 to m in increments of 1, where 
typically m=255. Stated differently [R,G,B] = [i,j,k] where 
i, j, k= 0-hti, where typically m=255. The output colors are 
produced according to output values produced from the input 
device stimulus data. A method for characterizing a color 
imaging system using the system 10 will be described in 
terms of the functionality of system 10. 

The processor 12 executes a software application 
program configured to carry out the method for 
characterizing a color imaging system. The software can be, 
e.g., in the C++ programming language. The processor 12 can 
be realized, for example, by a personal computer such as an 
Apple Macintosh™ or an IBM PC, or by a computer 
workstation. Alternatively, processor 12 could be realized 
by a microprocessor that accesses a read-only memory (ROM) 
into which the application program is loaded. The 
application program could be embedded in a color management 
software package, such as that provided with the Imation 
Corporation Rainbow™ color proofing system, commercially 
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available from Imation Corp., of Oakdale, Minnesota. The 
implementation of system 10 and the method using software 
affords flexibility in development and modification. The 
system 10 and method alternatively can be implemented by an 
integrated logic circuit for increased color processing 
speed. 

The color detector 14 and the color data memory 16 
allow the system 10 to obtain color data representing the 
output of various color imaging systems. The color data can 
be obtained directly from a color imaging system via the 
coior detector 14, or by accessing a color data file stored 
in the color data memory 16 . The color detector 14 and the 
memory 16 can help characterize, e.g., color plrinting 
systems, color display systems, color projection systems, 
and color-measuring systems. For use with a printing 
system, for example, color detector 14 is configured to 
measure color values of color patches 18 formed on the 
printing substrate 20. For use with a color display, color 
detector 14 is configured to measure color values produced 
on a phosphor screen or liquid crystal matrix. For use with 
a color projection system, color detector 14 is equipped to 
measure color values of projected colors. For use with a 
color-measuring system, color detector 14 is configured to 
measure color valves measured by the color-measuring system. 
As one illustration, color detector 14 could be configured 
to capture a scene or an animated sequence from a display or 
projection system, and generate color values representative 
of the captured imagery. 

The implementation of the color detector 14 depends 
on the application. The color detector 14 can be, for 
example, a color measurement system such as a Gretag™ SPM 
50 color measurement device, commercially available from 
Gretag, Inc., of Regensdorf, Switzerland, or a densitometer, 
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such as an X-Rite color densitometer, commercially available 
from X-Rite, of Grandville, Michigan. For display or 
projection system applications, color detector 14 can be a 
video camera or digital camera. 

The color detector 14 can store detected data in the 
color data memory 16 or in the processor 12. The color data 
obtained by color detector 14 can be loaded into color data 
memory 16 as a color data file, as indicated by line 22. 
Alternatively, the color data obtained by color detector 14 
can be directly loaded into a memory (not shown) associated 
with processor 12, as indicated by line 24. 

The processor 12 can access its own memory or the 
color data file stored in color data memory 16, as indicated 
by line 26, to obtain original color data previously 
detected by color detector 14. The color data memory 16 can 
store several color data files for a variety of different 
color imaging systems. Thus, a system user can direct 
processor 12, via a user interface 27 (e.g., a keyboard, a 
mouse, or a touch-sensitive screen) , to select one of the 
various color data files for purposes of characterizing a 
particular color imaging system of interest to the system 
user. 

The color data obtained from either color detector 
14 or color data memory 16 preferably represents CLE XYZ 
tristimulus values for each of a variety of color outputs 
generated by a color imaging system under study. 
Alternatively, the color data can be converted to CIE XYZ 
tristimulus values. The CIE XYZ tristimulus values 
represent the relative amounts of primary color stimuli 
required to match colors within a CIE color system. The 
relative values of XYZ are influenced by the power 
distribution of the illuminant, e.g., D 50 , and the CIE 
Standard Observer function, e.g., 2° or 10°. As an 
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alternative to CIE XYZ, the color data can be, for example, 
in the form of RGB data, CMYK density data, or other device 
dependent data. In addition, the color data preferably 
represents color outputs distributed across the color gamut 
of the color imaging system to provide a broad sampling for 
purposes of color characterization. Thus, in the case of a 
printing system, color patches 18 represent various 
combinations and gradations of colors formed by dye transfer 
or ink deposition. 

It was found that as colors, such as those of 
patches 18, on a source device varied from light to dark, 
and from neutral to saturated, the mismatch of colors 
produced on a destination device intended to match the 
source device colors also varied. Thus, even if light or 
neutral source device colors were reproduced on the 
destination device that matched well with the source device 
colors, dark or saturated colors did not match well, and 
vice versa. It was determined that in order to improve the 
match of both light and dark colors, and neutral and 
saturated colors, the color appearance model (s) used to 
convert the source device colors to destination device 
colors should vary with intensity and neutrality. The 
invention provides such a varying color appearance model . 

To characterize the color imaging system, processor 
12 converts the original color data using a color appearance 
model that uses a Y-reference and a white reference vector 
that depend on the intensities and neutralities of the 
original color data. In particular, the processor 12 varies 
the conversion of a Y value of the original color data to a 
luminance descriptor according the Y-reference value. The 
Y-reference value depends on the Y value of the original 
color data. The processor 12 varies the conversion of the 
original color data to two chrominance descriptors by 
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determining a white reference vector according to the 
neutralities of colors indicated by the trist imulus values. 
The term "white reference vector" refers to a vector of 
tristimulus values X r , Y r/ Z r that describe a color used as a 
white reference. 

The values of the Y- reference and the white 
reference vector are determined according to values of a 
"local white point" C LW and a "common white point" C^. The 
local white point C LW and the common white point are both 
vectors, with C LW = (X LW , Y LW , Z LW ) and C CT = (X^, Y^, Z w ) . As 
shown in FIG. 2, for determining the Y-reference and the 
white reference vector, values of the local white point C LW 
are used when the original color data indicates a near- 
neutral color and values of the common white point are 
used when the original color data indicates a non-neutral, 
i.e., saturated, color. Line 29, where proportions of X norm , 
Y norm , and Z^^ are equal, indicates neutral colors. The 
quantities X norra , Y norm , and Z norTn are the original color data 
normalized by reference white values such as from a 
perfectly white diffuser under standard viewing conditions. 
As shown by spherical surface 31 and tube 3 3 for purposes of 
illustration, there is an abrupt change between use of the 
local white point C LW and the common white point C LW . 
Preferably, the transition between use of the two white 
points is gradual, and not as abrupt as shown. 

The values of the vectors C LW and C^, depend on 
whether the device under consideration is absorptive (e.g., 
paper) or emissive (e.g., a monitor), and whether relative 
colorimetric intent or absolute colorimetric intent is 
selected, e.g., through the user interface 27. Depending on 
these factors, the local white point C LW and the common 
white point are equal to "media white," or "illuminant 
white, " or a predetermined white that may be independent of 
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the color imaging system under consideration. Media white 
for an absorptive medium is the white point of the 
absorptive substrate under a chosen illuminant. For an 
emissive device, media white is the whitest color that the 
device produces. Illuminant white is the color reflected by 
a perfectly diffuse white reflector illuminated by a 
selected illuminant . 

TABLES 1-4 illustrate the values of the local white 
point C LW and the common white point C cw for relative and 
absolute colorimetric intent for producing images using 
various combinations of absorptive and emissive media. As 
shown in TABLE 2, the local white point C LW for absolute 
colorimetry for two absorptive devices using different 
illuminants is the illuminant white corresponding to the 
illuminant used with each device. Also, the common white 
point for the case shown in TABLE 2 is the illuminant 
white of either medium, or another illuminant white, as long 
as the same white point is used for both media. As shown in 
TABLE 4, when two emissive devices are used the local white 
point for absolute colorimetry is the media white for each 
device. For two emissive devices the common white point 
shown in TABLE 4 can be a predetermined illuminant white, 
which may correspond to the media white of either emissive 
device. 

The local white point C LW for an emissive device 
using absolute colorimetric intent is the media white of the 
emissive device. Thus, as shown in TABLE 4, absolute 
colorimetric intent using two emissive devices merges with 
relative colorimetric intent because the media white point 
is used as the local white point for both rendering intents. 
It was discovered that when attempting to reproduce a color 
from an absorptive device on an emissive device using 
absolute colorimetry, the monitor produced a dark, overly- 
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yellow color. It was found that the human eye adjusts to 
define the media white of an emissive device as being white. 
By using the media white of the emissive device for the 
local white point C LW for absolute colorimetry, it was found 
that the emissive device can provide a good reproduction of 
the media white of the absorptive device and can provide 
good grey balance. The invention accounts for this finding 
by using the media white point for the local white point C fJM 
for absolute cclorimetric intent using emissive devices, 
thereby improving color match compared to prior techniques. 

Processor 12 is adapted to alter the color 
appearance model dependent on the neutralities and 
intensities indicated by the original color data to produce 
a modified L*a*b* color space, an I/a*b + color space. 
Processor 12 can store the I/a + b* color data in an I/a*b + 
memory 28, as indicated by line 30, in association with 
corresponding device stimulus data D. The LTa + b* color data 
based on color intensity and neutrality allow system 10 to 
achieve accurate color characterization over the entire 
range of a color imaging system. In particular, these data 
help ensure that system 10 can achieve substantial 
uniformity of the color characterization between light color 
shades and more intense colors, and between neutral and non- 
neutral colors. 

Referring to FIGS . 1 and 3, the functionality of 
system 10 is shown. The processor 12 can convert the 
original color data obtained from either color detector 14 
or color data memory 16 to converted color data. To, do so, 
the processor 12 uses any of the existing color appearance 
models, such as CIELAB, CIECAM-97s, RLAB, etc., subject to 
determination dependent on the intensities and neutralities 
of colors indicated by the original color data. The CIELAB 
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color appearance model will be used in the following 
description for purposes of illustration. 

Color data are obtained for a particular color 
imaging system by either detecting color data produced by 
the color imaging system, as indicated by block 40, or 
accessing color data in a color data file, as indicated by 
block 42. The color data have corresponding device stimulus 
data D used to actuate the color imaging system to produce a 
color corresponding to the color data for color-producing 
systems, and determined by measuring colors for color- 
measuring systems. The data indicated by block 44 comprise 
an array of device stimulus data D 0 -D m and CIE XYZ data t(X 0/ 
Y 0 , Z 0 ) . . . (X^, Y m/ Z m ) ] . 

The local white point C LW and the common white point 
are also obtained. As indicated by TABLES 1-4, 
depending on the color imaging devices involved and the 
colorimetric intent desired, one or more media white points 
and one or more illuminant white points are needed. The 
media white point (s) is (are) obtained from the data in block 
44. The illuminant white reference vector(s) is(are) 
obtained either by accessing a data file, as indicated by 
block 45, or by being provided as indicated by block 46 from 
an appropriate device, e.g., the processor 12, the color 
detector 14, or the user interface 27. These illuminant 
white reference vectors correspond to the illuminant white 
points indicated in TABLES 1-4 and discussed above. These 
vectors correspond to the , illuminant white of, e.g., the 
illuminant under which colors were viewed to produce the 
data array in block 44, or another similar block of data, or 
a predetermined illuminant white point. 

The color data indicated by block 44 are used to 
determine a white reference vector (X r ' , Y r ' , Z r ' ) and a Y- 
reference Y r " in relation to the local white point and the 
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common white point, as indicated by block 48. The white 
reference vector (X r ' , Y,' , Z r ' ) and the Y-reference Y c " are 
determined for each set of color data shown in block 44 
according to the intensity and neutrality of the color 
indicated by the color data. The white reference vector 
(X r ' , Y r ' , Z r ') and Y-reference Y r " are determined according 



to the following equations: 

Y r " = Y.j, (1 - sat(Y,Y._ w ) ) + Y CT * sat (Y,Y LH ) [6] 
sat ( Y , Y LW ) - 1.0 - (Y/Y LW ) [7] 
X r ' = X La (1 - sat(C,C LW )) + Xa, * sat(C,C LW )) [8] 
Y r « = Y La (1 - sat(C,C w )) + Yc * sat (C, C w ) ) [9] 

Z r ' = Zu, (1 - sat (C,C,_.J ) + Z CT * sat(C,C LW )) [10] 

sat(C,C LW ) = (devX' Y' Z ' /maxDev) r [11] 

maxDev = sqrt (6 . 0/9 . 0) x max ( X ' , Y ' , Z ' ) [12] 
devX'Y'Z' = sqrt((X'- avgX'Y'Z') 2 + (Y' -avgX' Y' Z' ) 2 

+ (Z' -avgX'Y'Z' ) 2 ) [13] 

avgX'Y'Z' = (X' + Y' + Z')/3.0 [14] 



where X'=X/X LH/ Y'=Y/Y LVI , and Z'=Z/Z L „, the color vector is 
C=(X,Y,Z) with X, Y, and Z being a set of color data from 
block 44. As indicated by block 50, the Y-reference Y r " and 
the white reference vector (X r ' , Y r ' , Z r ' ) are output for 
further use. 

As indicated by block 52, the white reference vector 
(X/ , Y r ' , Z r ' ) and Y-reference Y r " are used to determine the 
L*a*b* color space by converting the original color data from 
block 44 according to the following color appearance model 



equations : 

1/ = 116 x f (Y/Y r " ) - 16 [15] 

a* a 500 x [f(X/X r ') - f(Y/Y r ')] [16] 

b* = 200 x [f(Y/Y r ') - f(Z/Z r ')] [17] 

f (cj) = (w) l/:i (J>0. 008856 [4] 
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f (o>) =7 .787 (a>) +16/116 ojsO . 0 08 856 



[5] 



Applying equations [15] -[171 to the color data 
obtained for the color imaging system, in combination with 
determining the white reference vector (X r ' , Y r ' , Z r ' ) and Y- 
reference Y r ,! as in equations [6] - [14] , produces a set of 
LTa^b* color space data that characterizes the color response 
of the color imaging system. The L/a + b + color space data are 
stored in an array as indicated by block 54 . 

The lightness value L* depends on the Y values of 
the color data in block 44 and of the local white point and 
common white point. Equations [6] - [7] and [15] indicate 
that as the value of Y in the original color data from block 
44 increases, corresponding to a lighter color, the Y- 
reference Y r " approaches the Y value for local white. As 
the value of Y decreases, corresponding to a darker color, 
the Y-reference Y r " approaches the Y value for common white. 
Thus, the higher the Y value in the set of original color 
data from block 44, the more the lightness value L* depends 
on the local white point as a reference white point and the 
lower the Y value, the more the 1/ value depends on the 
common white point as a reference white point. 

The redness-greenness value a + and the yellowness- 
blueness value b* depend on the X, Y, and Z values of the 
original color data in block 44 and of the local white point 
and common white point. Equations [8] - [14] and [16] - 
[17] indicate that as the color indicated by a set of 
original color data X,Y,Z, approaches a neutral color, the 
white reference vector (X r f , Y r ' , Z r ' ) approaches the local 
white point C LW . As the color indicated by the set of color 
data deviates from a neutral color, the white reference 
vector (X r # , Y r ' , Z r ' ) approaches the common white point C^. 
Thus, the more neutral the color indicated by the set of 
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original color data from block 44, the more the redness- 
greenness value a' and the yellowness -blueness b + value 
depend on the local white point as a reference white point 
and the more saturated the color, the more the a' and b + 
values depend on the common white point as a reference white 
point . 

As shown in FIG. 4, a system 32 for performing a 
color transformation includes a processor 12, a color 
detector 14, a color data memory 16, a first substrate 20, a 
second substrate 38, and a stored color transformation table 
34. The components of system 32 substantially correspond to 
system 10 of FIG. 1, with similar features having identical 
reference numbers. System 32 performs similar functions to 
those of system 10 to characterize two or more different 
color imaging systems. System 32 calculates a mapping 
between the color imaging systems based on the 
characterizations. The mapping is used to generate the 
color transformation table 34, which is stored in memory by 
the processor 12 as indicated by line 36. The color 
transformation table 34 is used to produce a color on a 
target color imaging system that visually matches the color 
produced by a source color imaging system. 

The color detector 14 is configured to obtain color 
data representing the output of two or more color imaging 
systems. As shown in FIG. 4, color detector 14 is 
configured to detect color data from color patches 18 formed 
on the first printing substrate 20 by, or provided to, a 
first color imaging system and color patches (not shown) 
formed on the second printing substrate 3 8 by, or provided 
to, a second color imaging system. Alternatively, detector 
14 can be configured to measure color values produced by 
various emissive or projecting devices, such as by phosphor 
screens or liquid crystal matrices. The detector 14 can 
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send the detected color data to the color data memory 16 
along line 22 for storage or to the processor 12 along line 
24 for storage or processing. 

The color data memory 16 stores color data in one or 
more color data files. The color data can be received from 
the detector 14 or from another source. The color data 
memory 16 can transmit stored color data co the processor 12 
along line 2 6 . 

The processor 12 is configured to receive color data 
from the detector 14 and/or the memory 16 and to process the 
received data to perform a color transformation. The color 
data detected by the detector 14 and/or stored in color data 
memory 16 preferably represents CIE XYZ tristimulus values 
for each of a variety of color outputs generated by the 
different color imaging systems. To perform a color 
transformation, processor 12 can convert the color data 
obtained for each of the different color imaging systems 
using the I/a'b + color appearance model. The processor 12 
can convert original color data obtained for the first and 
second color imaging systems into L*a + b* data. The 
conversion uses white reference vectors and Y-reference 
values that vary between local white and common white values 
of each color imaging system. The Y-reference values and 
white reference vectors depend on the intensities and 
neutralities of the colors indicated by the original color 
data from each color imaging system. 

FIG. 5 illustrates the functionality of system 32 
and a method for performing a color transformation. Color 
data are obtained for a first, or source, color imaging 
system, as indicated by block 56, and for a second, or 
target, color imaging system, as indicated by block 58. The 
resulting color data for the first color imaging system 
comprises an array of CIE XYZ data [(X 01 , Y 0#1 , Z 0<1 ) . . . 
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(X^!, Y ml , Z ml )], whereas the color data for the second 
color imaging system comprises an array of CIE XYZ color 
data t(X 0f2 , Y 02 , Z 0 2 ) . . . (X m 2/ Y m>2 , Z fflj2 ) ] . For example, 
the vector (X 0 : , Y 01 , Z 0>x ) represents the color data 
obtained for the media white point associated with the 
first, absorptive (e.g., CMYK), color imaging system, and 
the vector (X 0 2 , Y 0 2 , Z 0 2 ) represents the color data 
obtained for the media white point associated with the, 
absorptive, second color imaging system. Alternatively, 
these vectors can represent the media black points for 
emissive systems such as RGB systems. Similarly, the vector 
(X m fl , Y m#xi Z m J represents the color data obtained for a 
maximum (minimum) intensity color produced by the first 
absorptive (emissive) color imaging system, and the vector 
(^.2/ Y m.2' z m,2) represents the color data obtained for a 
maximum (minimum) intensity color produced by the second 
absorptive (emissive) color imaging system. The value of 
X nt/1 is given by the value of X x which satisfies: 



MAX(ABS(X 1 - X b ) ) 

for all values of X x measured from all possible colors 
imaged by the 1 first color imaging system, where ABS returns 
the absolute value of the difference value in parentheses, 
and MAX returns the maximum difference value for all values 
of X x relative to X b ; Values of Y m#1 , Z m l , X^, Y m , 2 , and Z TOf2 
are similarly defined. 

Local white points C LW1 and C LW2 and the common white 
point Co, are obtained from the appropriate media and 
illuminant white points (TABLES 1-4) . The media white 
points are obtained from the data in blocks 56 and 58 and 
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the illuminant white point (s) is (are) obtained by accessing 
stored data (block 59) or by being provided (block 60) as 
appropriate. The illuminant white points correspond to, 
e.g., the illuminant under which colors were viewed to 
produce the data in blocks 56 and 58, or a predetermined 
white point. 

As indicated by block 62, the color data in blocks 
56 and 58 are used to determine white reference vectors 
(X rx ', Y rl ', Z rl ') and (X r2 \ Y r2 ' , Z r2 ') and Y-stimuli Y rl " and 
Y r2 " in relation to local white points C LW1 and C LW2 and the 
common white point C^. These values are determined 
according to intensities and neutralities of colors 
indicated by the data in blocks 56 and 58 using the 



following equations: 

Y rl M = Y LW1 (1 - sat ( Y , Y LW1 ) ) + Y^ * sat (Y,Y LW1 ) [18] 

sat(Y,Y LW1 ) - 1.0 - (Y/Y LH1 ) [19] 

X rl ' = X LW1 (1 - sat (C 1# C LW1 ) ) + X CT * sat(C lf C LW1 ) ) [20] 

Y rl ' = Y LW1 (1 - sat(C 1/ C LW1 ) ) + Ye * sat (C, , C LW1 ) ) [21] 

Zn' = Zlhi (1 - sattC^C^)) + * sat (C 1# C LW1 ) ) [22] 

sat (C lt Cuti) (devX 1 'Y 1 'Z 1 '/maxDev) Y [23] 

maxDev = sqrt (6 . 0/9 . 0) x max(X x ' , Y x ' , Z x ' ) [24] 
devX/Y/Z/ = sqrt((X 1 # - avgX, ' Y x 1 Z l ' ) 2 + (Y/ -avgX L ' Y/ Z/ ) 2 

+ (ZZ-avgX.'Y/Z/) 2 ) [25] 

avgX/Y/Z,' = (X/ + Y x ' + Z L ')/3,0 [26] 

Y r2 " = Y LW2 (1 - sat(Y,Y LW2 )) + Y^ * sat (Y,Y LH2 ) [27] 

sat(Y,Y LW2 ) = 1.0 - (Y/Y LW2 ) [28] 

X r2 ' = X LW2 (1 - sat (C 2 ,C LW2 ) ) + X OT * sat(C 2/ C LW2 )) [29] 

Y r2 ' - Y LW2 (1 - sat (C 2 , C LW2 ) ) + Yew * sat (C 2 ,C LH2 ) ) [30] 

Z r2 ' = Z LH2 (1 - sat(C 2 ,C LH2 ) ) + Z OT * sat (C 2 ,C LW2 ) ) [31] 

sat (C 2> C LW2 ) = (devX 2 'Y 2 'Z 2 '/maxDev) Y [32] 
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maxDev = sqrt (6 . 0/9 . 0) x max(X 2 ' , Y 2 ' , Z 2 ' ) [33] 
devX 2 'Y 2 'Z 2 ' = sqrt((X 2 '- avgX 2 ' Y 2 ' Z 2 ' ) 2 + (Y 2 ' -avgX 2 ' Y 2 ' Z 2 ' ) 2 
. + (Z 2 ' -avgX 2 'Y 2 'Z 2 ' ) 2 ) [34] 
avgX 2 'Y 2 'Z 2 ' = (X 2 ' + Y 2 ' + Z 2 ')/3.0 [35] 



where X x ' ^X x /X L]filt Y x ' =Y 1 /Y LW1 , and Z 1 '=Z 1 /Z LW1/ the color vector 
C t = (X lf Y lf Z x ) with X 1# Y 1# and Z, being a set of color data 
from block 56, and where X 2 '=X 2 /X LW2 , Y 2 '=Y 2 /Y LW2 , and 
Z 2 , =Z 2 /Z LW2/ the color vector C 2 = (X 2 , Y 2# Z 2 ) with X 2 , Y 2/ and Z 2 
being a set of color data from block 58. As indicated by- 
blocks 64 and 66, the values for Y-stimuli Y rl " and Y r2 " and 
for white reference vectors (X rl ' , Y rl ' , Z rl ') and (X r2 ' , Y r2 ' , 
Z r2 ' ) are output for further use. 

As indicated by block 68, the white reference 
vectors (X rl ' , Y rl ' , Z rl ') and (X r2 ' , Y r2 ' , Z r2 ') and Y-stimuli 
Y rL " and Y r2 " are used to determine the L^a^b^ and L 2 *a 2 + b 2 + 
color spaces by converting the original color data from 
blocks 56 and 58 according to the following color appearance 



model equations:. 

V = US x f^/Y^") - 16 [36] 

a x -= 500 x [f(X l /X rl ') - f(Y x /Y rl ')l [37] 

b x * = 200 x [f(Y x /Y rl ') - f(Z x /Z rl ')] [38T 

L 2 * = 116 x f(Y 2 /Y r2 ") - 16 [39] 

a 2 * = 500 x [f(X 2 /X r2 ') - f(Y 2 /Y r2 ')] [40] 

b 2 + = 200 x [f (Y 2 /Y r2 ' ) - f(Z 2 /Z r2 ')] [41] 

f (w) =(cj) 1/3 w>0. 008856 [4] 

f (w) =7 .787 (cj) +16/116 cjsO. 0088 56 [5] 



Applying equations [36] - [41] to the color data 
obtained for the color imaging systems, in combination with 
determining the white reference vectors (X r / , Y rl ' , Z rl ') and 
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(X r2 '/ Y r2 ', Z r2 ') and Y-stimuli Y rl " and Y r2 " as in equations 
[18] -[35], produces sets of L 1 *a l *b l * and L 2 *a 2 *b 2 * color space 
data. These sets of data are subjected to mult i -dimensional 
interpolation by nhe processor 12 to completely characterize 
the color responses of the color imaging systems. The 
complete sets of L 1 *a 1 *b x * and L/a/b 2 * color space data are 
stored in arrays as indicated by blocks 70 and 72 
respectively. 

The complete sets of En/a^b^ and L 2 *a 2 *b 2 * color space 
data are mapped to each other by the processor 12 as 
indicated by block 74. The processor 12 maps the color 
space data from the two color imaging systems to each other. 
This also maps the device stimulus data corresponding to the 
color space data for each color imaging system to each 
other. Thus, a mapping between device stimulus data (e.g., 
CMYK to C'M'Y'K' or RGB to CMYK) for the two color imaging 
systems is produced. Using this mapping, once the stimulus 
data associated with one color imaging system when 
producing, or measuring as the case may be, a color are 
known, the other system can be actuated with the mapped 
stimulus data to reproduce the color. The mapping is stored 
in a color transformation table as indicated by block 76. 

Other embodiments are within the scope of the 
appended claims. For example, color data can be obtained in 
manners other than those described. Color appearance models 
that employ a luminance descriptor and descriptors for the 
relative amounts of red, green, yellow, and blue, other than 
those mentioned, can be used. Also, although the XYZ 
tristimulus values were described above as having, e.g. ,. 256 
discrete values (e.g., for a digital system), they may take 
on values for an analog system over a broad range of values 
to characterize, e.g., paint samples. 
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It has been found that calculating the white 
reference vector for dark colors based upon the angle from 
neutral may be helpful. Dark colors such as deep blue may 
have large chromas (e.g., about 60.0) but low values of XYZ 
and therefore devXYZ . This would result in the dark 
saturated colors being referenced to local white (or media 
white) rather than common white. Referencing these colors 
to common white, however, yields better color reproduction. 
It has been found to be helpful in accurately characterizing 
and reproducing colors to use a value for deviation from 
neutral that is based on the angle from neutral or on a sine 
or tangent of that angle, although the invention is not 
limited to these values of deviation from neutral. Appendix 
A provides an example, which in no way limits the scope of 
the invention or the appended claims, of C++ code utilizing 
the tangent of the angle from neutral. The code in Appendix 
A could be modified to use the angle, or the sine of the 
angle from neutral . 

Other relationships, than provided above, between 
I/a*b + values and the local white point C LW and the common 
white point C cw are possible. The relationship between the 
white reference vector (X r ' , Y r ' , Z r ' ) and the Y-reference 
and the local white point C LW and the common white point 
was exemplary only and other relationships are possible. 
For example, the L*a*b* color space can be determined 
according to the following equations. 

L* = (1.0 - sat_L*) x L* rel + sat_L* x L* ab3 [42] 
(a*,b + ) = (1.0 - sat_a*b*) x (a*,b*) rel + sat_a*b* x 

(a*,b*) ab3 [43] 

where 

sat_L* - 1.0 - ( Y/Y LW ) [44] 
sat_a*b* = C*/L* [45] 
C = sqrt(a* 2 + b* 2 ) [46] 
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Here the scaling from local white to common white is 
performed by linear addition of conventional absolute and 
relative CIELA3 calculations given by equations [i] - [5] . 
Other relationships, such as nonlinear scaling, are 
possible. 

Also, the Bradford chromatic adaptation transform 
can be modified to scale white points when using absolute 
white points instead of normalized white points such as for 
media-relative mapping of reflective prints. It is useful 
to think in terms of reference conditions (the destination 
conditions) and local conditions (the source conditions) . 
For example, in the case of a reflection print, the white 
point for the reference conditions is usually the reference 
illuminant (e.g., D50) , and the white point for the local 
condition is the media white. 

First, the XYZ values are normalized using the 
luminance (Y) and the local black point [X 1)c/ Y lk , Z lk ] and 
then transformed through the Bradford matrix M b to obtain 
modified responses of the cones of the human eye: 
TrJ r(X 1 -X lk )/(Y 1 -Y lk )1 
|Gj- Mj (Yi-Y lk ) / (Yi-Y lk ) | 
LBj L(Zi-Z r J/(Y 1 -Y lk )J 

[ 0.8951 0.2664 -0.16141 
M b =| -0.7502 1.713 5 0.03 67| 
L 0.0389 -0.0685 1.0296J 

where X 1# Y x , and Z, are the tristimulus values of a color to 
be transformed, and R l9 G 1# and B x are the modified responses 
of the cones of the human eye. 
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Second, the RGB signals are chromatically adapted 
from the local condition to the reference condition using 
the following equations (assuming complete adaptation) : 



R re f — (Rrv/R-iw) x R-l 

G ref = (G„/G lw J x G x 

B ref = Sign[B l ]x(B r ,/B 1 ,/)x|B 1 |* 

]8 = (B U /B W ) MIM 

where R r „, G r ,, and B rrf are the reference white RGB values, 
R lw , G lw , and B Xw are the local white RGB values. 

Third, the adapted RGB signals are transformed back 



where is the reference white luminance value. 

This adaptation to the reference condition takes 
into account any differences between the absolute level of 
the reference and local white luminance, in addition to 
factoring in the local black point. By scaling the Y values 
based on the difference between local and reference 
conditions, this adaptation method can be used for media- 
relative proofing to map between reflection prints. In this 
application, it is undesirable to simulate the absolute 
white point by printing color in the background area of the 
proof. 



(neutral) axis is mapped to the neutral axis of the 
reference condition, while larger adaptations are applied to 



to XYZ: 



rx ref i 

|Y reC | = M b - 1 

fZr.fl 



FR ref x Y x x Y rw /(Y lw - Y lk )l 
|G W . x Y l x Y w /(Y lw - Y lk ) I 
[B ref x Y x x Y rw /(Y lw - Y U )J 



In addition, the local white to local black 
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colors that are further away from the local neutral axis. 
In this context, this can be viewed as an alternative 
approach to the variable white reference concepts discussed 
above . 

The XYZ values that have been adapted to the 
reference conditions are used to compute L'a'b* values: 

LT = 116 x f(Y/Y n ) - 16 

a* = 500 x [f (X/X n ) -f (Y/YJ ] 

b' = ,200 x [f (Y/Y n ) -f (Z/Z n ) ] 

f (u) = (co) 1/3 u>0. 008856 

f (cj) =7 .787 (w) +16/116 cjsO. 008 856 

This chromatic adaptation formulation can use color 
spaces other than the standard I/a*b fc color space. This 
formulation can just as easily be used with alternative 
color spaces such as the ones from CIECAM-97s, LLab, or 
RLab. By utilizing the same XYZ to "uniform color space" 
transformation for all components of a given matching step, 
however, different adaptation methods can be applied to 
achieve useful results, as discussed below. 

The concepts described can be applied to chromatic 
adaptation methods other than the Bradford adaptation, such 
as the von Kries transform. Thus: 

rx ref l [L^ 0 0l oi rxj 

|Y ref |= M/ 1 I 0 M„ 0|| 0 l/(M lw -M lk ) OlMjYj 

LZreJ L 0 0 S J L 0 0 1/ ( J LzJ 

r 0.38791 0.68898 -0.07868] 
M^I-0. 22981 1.18340 0.0464l| 
L 0 0 1.0 J 
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where [L r ,, M rw , S^J are the LMS (long, medium, and short 
wavelength band) values of the reference white, [L lw , M lw , 
3 lw ] are the LMS values for local white, and [L 1)c , M lk , S lk ] 
are the LMS values for the local black. The resulting 
adapted XYZ value is used as the starting point for 
computing the L*a*b* value using the XYZ to Lab equations 
described above. 

These modified Bradford and von Kries adaptations 
can be applied differently than conventionally thought. 
These adaptations can be used in conjunction with 
standardized CIE-Lab (e.g., Hunt, LLab, RLab, CIE-CAM97) 
calculations and can be applied to different types of media 
(e.g., CRT, reflection print) when performing cross-media 
rendering. To perform such rendering, colors from each 
media are converted using a modified transform, including 
reference and local XYZ and/or LMS values, and a desired 
color space so that the converted colors can be mapped. The 
following are three specific examples of ways in which these 
methods can be applied: 

1. Matching print to print using media-relative 
colorimetry : 

XYZ^ is the XYZ of the reference illuminant for both 

print devices, assuming a common illuminant. 

XYZ lw for each of the prints is the actual XYZ of the 

media white of the respective print. 

Bradford-type adaptation is used for both prints. 

In this case, the illuminant white will be common 
for the two prints and will reflect the desired viewing 
environment (e.g., a D50 light booth). White on the first 
print will map to white on the second print, and differences 
in the background paper color will be taken into account in 
generating the match for chromatic colors. 
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2. Matching print to display using media-relative 
colorimetry : 

XYZ r „ is the XYZ of the reference illuminant. 
XYZ lw for the print is the actual XYZ of the media 
white of the print. 

XYZ,., for the display is the actual XYZ of the 
monitor white point of the display. 

Bradford- type adaptation is used for both the print 
and the monitor. 

In this case, the monitor colors will be 
chromatically adapted to the reference illuminant so as to 
give a good perceptual match. By scaling the luminance in 
the Bradford equations to account for differences between 
the Y of the media white and the Y of the illuminant white, 
the white point of the print will map to the white point of 
the monitor. 

3, Matching print to display using absolute 
colorimetry: 

XYZ rw is the XYZ of the reference illuminant. 
XYZ lw for the print is the actual XYZ of the media 
white of the print. 

XYZ lw for the display is the actual XYZ of the 

monitor white point of the display, 

Bradford-type adaptation is used for the monitor. 

Absolute CIE-Lab is used for the print. 

By performing matching in this way, instead of the 
standard approach to this situation of using absolute 
colorimetry for both the monitor and the print, the 
Bradford-type adaptation performed on the monitor values 
simulates the chromatic adaptation occurring in the eye. 
Absolute colorimetry, performing no chromatic adaptation, 
can be used on the print while performing the chromatic 
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adaptation on the monitor to achieve a good perceptual 
match. This approach could be applied using a different 
conversion between XYZ values and the matching color metric 
(e.g., L.a.b. ( LLab, RLab, CIECAM-97s) . 

While the descriptions of these three exemplary 
applications assume no black point mapping, so XYZ lk = 
(0,0,0), for all of these applications media black point 
mapping can be applied. Media black point mapping is 
described in patent application serial number 08/884,411 
filed June 27, 1997, entitled "Characterization of Color 
Imaging Systems," assigned to the assignee of this 
application, and incorporated herein by reference. 

What is claimed is: 
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